Phase separation in double exchange systems 
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Ferromagnetic systems described by tlie double exchange model are investigated. At temperatures 
close to the Curie temperature, and for a wide range of doping levels, the system is unstable toward 
phase separation. The chemical potential decreases upon increasing doping, due to the significant 
dependence of the bandwidth on the number of carriers. The reduction of the electronic bandwidth 
by spin disorder leads to an enormously enhanced thermopower which peaks near Tc, with a sign 
opposite that predicted by a rigid band model. 
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Doped manganese oxides exhibit many unusual fea- 
tures, most notably the eponymous phenomenon of colos- 
sal magnetoresistance (CMR) jl]-^, and a phase transi- 
tion from a high temperature paramagnetic insulator to a 
low temperature ferromagnetic metal. The Mn 3d states 
split into a lower t2g triplet, forming an 5* = | core spin, 
and an upper eg doublet, the conduction band. This 
physics is described by the double exchange model [Q: 
the conduction electrons which hop throughout the lat- 
tice are ferromagnetically coupled to the local core spins 
because of Hund's rules. Typical values for the Cg hy- 
bridization t and intra-atomic exchange Jh are t ^ 0.1 eV 
and Jh ~ 1 — 3 eV. Additionally, one may choose to incor- 
porate other terms, most notably Heisenberg couplings 
between neighboring core spins (due to superexchange) 
and Jahn- Teller (JT) phonons which break the degener- 
acy of the Cg level. 

The simplest model, however, is one in which the de- 
generacy of the Cg orbital is simply ignored (possibly due 
to a cooperative JT distortion), and Jh is set to infinity. 
In this limit, the spin of each eg electron must agree with 



that of its t2g core, and we may write e 



t - J 



where 



Zjg. is the spinor describing the orientation of the core 
at site i: z^^ — cos(^^i), z.^ ~ sin(i0j) exp(— z0j). The 
Hamiltonian is then 



H = -t 



H.c. 



(1) 



The hopping is maximized when neighboring core spins 
are parallel. The ferromagnetic interaction induced in 
this way describes qualitatively the physics of the man- 
ganites. The Curie temperature increases with the num- 
ber of carriers, and, in the paramagnetic phase, the ten- 
dency towards localization is enhanced. On the other 
hand, it has been extensively argued that the double ex- 
change model does not suffice to describe the insulating 
behavior at high temperatures ||]. A disordered distri- 
bution of spin gives rise to off-diagonal disorder, which 



is in a sense weaker than the more conventional diag- 
onal disorder, and most of the states in the band are 
delocalized [|| (see, however, the comments of 0). The 
eg bandwidth remains finite and is unremarkable as one 
passes through the Curie temperature Tc Scatter- 
ing by magnetic fluctuations near Tc also is insufficient 
to explain the insulating regime A variety of phe- 
nomena, including variable range hopping Q, formation 
of spin polarons |jlffl , and the dynamic Jahn- Teller effect 
|]ll| have been invoked in an attempt to understand the 
insulating paramagnetic phase. 

Here we will show that the double exchange model (|l|) 
is unstable toward charge segregation at low to moder- 
ate doping and temperatures near Tc. This phenomenon 
arises due to the coupling between magnetic fluctuations 
and the electronic chemical potential. Exchange interac- 
tions between core spins tend to suppress this effect, but 
do not completely eliminate it. The Coulomb interaction 
will suppress complete phase separation of charge carri- 
ers, resulting in microscopic domains of different charge 
density. This inhomogeneity in turn enhances localiza- 
tion of electrons. A small ferromagnetic region of high 
density, surrounded by a low density paramagnetic back- 
ground, is indistinguishable from the spin polaron picture 
discussed in the literature. 

To see why phase separation should occur i n (pi ), let 
us reexamine de Gennes' mean field approach |12|| 7 The 

quantity W = ^/{}^~^~^ - (cos2(i^y))i/2 is a "spin 
reduction factor" which multiplies the fermion hopping 
t, compressing the bare {W = 1) dispersion E to WE. 
At zero temperature, states with E < Ep are filled, and 
fj. = WEp. The density of states D{E) and doping x = 
1 

2 



S determine fj. via <5 = ^dE D{E). Thus, 
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H dW 
W 'W 



(2) 



When dfi/dS < 0, the system can lower its energy 



1 



through phase separation. Now at zero temperature we 
should expect perfect order in the spins, hence W = 1 
and the compressibiUty is positive. At finite tempera- 
ture, however, we expect fidW/dS < since increased 
carrier density \6\ enhances the fermion kinetic energy, 
which serves as an exchange couphng for the spins (this 
is true whether the carriers are electrons or holes) . If we 
assume the fermions to be degenerate, the condition for 
negative compressibility is ^iD{p,/W) dW / dS < —W . 
To satisfy this relation requires both finite T and low to 
intermediate \6\ (note /i = for S — 0, assuming a sym- 
metric band). As we shall see, this behavior may per- 
sist above the Curie temperature as well. At sufficiently 
large T, however, the fermions are nondegenerate and 
the chemical potential dominates the bandwidth, giving 
(5 = i tanh(/i/2fcBT) and dn/dS = ^KT/{1 - 4(5^) > 0. 

This picture also leads to an anomalous thermopower 
psf . At low temperatures the variation of the chemical 
potential with T arises from two sources: the Fermi dis- 
tribution {e.g. the Sommerfeld expansion) and the band 
narrowing W . Thus one obtains 
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dW 



fJ-o 



dT 3 D(Aio) 



KT 



(3) 



For models where D'{E)/E > 0, such as the elliptical 
density of states, the first term, arising from the band 
narrowing, opposes the second. As we shall see for all but 
the lowest and highest temperatures it is the first term 
which dominates, and the critical behavior of W in the 
vicinity of leads to a peak in the thermopower which 
is several orders of magnitude larger than that usually 
encountered in metals. 

Adding purely ferromagnetic superexchange between 
the core spins reduces the dependence of on 5 and 
thereby opposes phase separation. In Lai_2:A2;Mn03, 
however, the interplane superexchange is antiferromag- 
netic for small x, and its competition with double ex- 
change leads to a canted structure ||l^. The bandwidth 
then depends on the canting angle and increases with in- 
creasing number of carriers. This property leads to phase 
separation even at T = [p^[-p^. 

There is a subtlety in the mean field theory surround- 
ing which density of states D{E) one should use. Inte- 
grating out the fermions from (pi) gives a free energy 



F 



E 

p=i 



Tt HP -TS[{fli\] (4) 



where /(e) = -fcBrin(l + e-"/^^^) and H,.j 
with Tij = t li i and j are nearest neighbors and other- 
wise. The entropy = —k^ T^r q\tiq and the trace 
in (^ are computed from a trial density matrix q for 
the spins. Assuming an uncorrelated g — P{fli) with 
P(!7) oc ex^iQ^z), we have that {ex-p{i(f>i)) — 0, hence 
in the locator expansion of Tr^y^ j H^, many of the 2^ 
terms associated with a given length p path average to 
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FIG. 1. Phase diagrams (a) for the mean field theory of 
eqns. (|,§), and (b) for that of eqns. (^,^), both computed 
using an elliptic density of states. The thick solid line is Qc{S) . 
The thick dashed line is Ss(5), the boundary of the region of 
phase separation. In (b), the thin solid line denotes 
the temperature where the bandwidth collapses. 



zero. These cancellations are avoided for paths which re- 
trace themselves, for which each z^^z^^ term has a z^^z^^ 
mate. Thus, in the vicinity of Tc where magnetic fluctu- 
ations are significant, it is better to use the retraced path 
approximation jT^ , 



i^(7) = -^ 



4 ? 
-.1 



(5) 



where z is the lattice coordination number, in one's mean 
field calculations. £'(7), given above with 7 = —E/B in 
units of the bare half-bandwidth B = 2{z — lY^^t, in- 
terpolates between the exact one-dimensional tight bind- 
ing density of states at 2; = 2 to a semi-elliptic form at 
z ~ 00. The hopping is then modulated by W, the root- 
mean-square average of Zj^^Zj^. 

The mean field equations of de Gennes' model are 



(6) 



where M — ctnh (Q) — Q^^ is the magnetization, and 
The temperature and chemical po- 



w 

tential are scaled by a = 11/ B and 8 = k^T/B. 

The mean field solution exhibits a Curie temperature 
&c{S). With t w 100 meV the cubic lattice model 
has B = 2^/5 1 = 5200 K. We find a maximum 9c = 
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0.0545 at (5 = 0, corresponding to a Curie temperature 
of Tc = 280 K, very much consistent with experimental 
values. M vanishes at 8c, and for 8 > 0c the bandwidth 
is finite, reduced by the spin factor W = The com- 
pressibility = da/ 86 and thermopower S = da/dQ 
are discontinuous at the transition. S peaks near Tc with 
an anomalously large value; we shall return to this point 
below. For the z = oo case, we find a crescent region 
of phase separation 9 e [6s(5), 0c(5)] extending from 
6 = 0.382 to (5 = i, in qualitative agreement with our 
earlier discussion. In Fig. ^(a) we plot the mean field 
phase diagram for the limiting case z — oo. (As one ap- 
proaches the one-dimensional limit z = 2, the region of 
phase separation is enlarged.) For z > 2 we find 0c oc x 
and Qs oc x^^^ asa; = ^ — ^^0. (The model is of course 
invariant under x x.) 

An unphysical aspect of the preceding mean field anal- 
ysis is that W, which should be determined by local spin 
correlations, is tied to the magnetization M. Initially W 
decreases linearly with 6, leading to artifactual behavior 
in the low temperature magnetization and thermopower. 
The root of the problem is the classical treatment of the 
core spins. At temperatures 9 < 0{1/S), one must treat 
the spins quantum mechanically. We should then expect, 
from spin wave theory, that W = 1 ~ 0(8^/^). 

A somewhat more sophisticated mean field theory 
can be implemented using the Schwinger boson p7| ] 
method. A nondynamical field enforces the constraint 
^itT-^io- ~ 1 at every site. The core spins are quan- 
tized according to [z^^^z^^,] = SijSaa' /^S . The mean 
field Hamiltonian is obtained through a Hartree decou- 
pling of the bosonic z^^z- and fermionic c^c - hopping 



terms: 

= N{ztWK - A) - ^ 



c^c, + A^; 



(8) 



where K — {c'-c^ 



^ " ^ and W — l^z^^z-^) Such a model was in- 
troduced by Sarker |20[ , who identified a Curie transition 
and found that the eg fermion band becomes incoherent 
above Tc. For our purposes we are interested in phase 
separation. Accounting for the possibility of condensa- 
tion of Schwinger bosons, we write = {z^,^)- Assum- 
ing condensation only at fc = 0, we define p = \^k=o,cy\^- 
The mean field equations are then 
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W = p + 



-y_d7 7i^(7)ctnh(^^ 

1 

(5 = - / djD{-f) tanh 
2 7-1 

1 



c?7 7 D(-/) tanh 
1 



a + 
29 

a + W"/ 
29 



(9) 



where A = X/zt, 9 = k^T/zt, and a = p/zt {B = zt). 

We have solved the mean field equations for a semi- 
elliptic density of states. Several features of the solution 
are noteworthy. Again as expected there is a Curie tran- 
sition at 9c ((5). In the absence of an external magnetic 
field, the magnetization M is equal to the condensate 
fraction p and vanishes for 9 > 9c ((5). Second, the quan- 
tity W ^ which is directly proportional to the c fermion 
bandwidth, is unremarkable through Tc. It decreases 
monotonically and eventually vanishes at a temperature 
9* (5) = (1 + 5-1)1/ Vl - 4^2/8^2. rpj^g dimensionless 
electronic kinetic energy K, which is proportional to the 
boson bandwidth, also vanishes at this point. The van- 
ishing of K and at 9* is likely a spurious artifact 
of the mean field theory; it was noted in ||l^,|l6|. How- 
ever, in a range 9c < 9 < 9* above the Curie transition 
this theory describes a state with vanishing magnetiza- 
tion yet finite and 9-dependent fermion bandwidth. The 
mean field parameters are shown versus 9 in Fig. |. 

Third, we find that k~^(Q) has a discontinuity in slope 
at 9c and a jump discontinuity at 9*. The locus of points 
where k~^{6,Q) — marks the boundary of the region 
of phase separation. This is shown in Fig. |](b). Again 
corroborating our initial discussion, we find that phase 
separation occurs for \5\ > Sq ~ 0.289, and for a range 
of temperatures 9 S [9s((5), 9*((5)] surrounding 9c. The 
upper boundary is 9* ((5) itself; this is likely an artifact of 
the mean field theory. We expect that W should continue 
to decrease as 9 increases, tending to a finite asymptotic 
value that this mean field approach cannot describe. At 
any rate, for large 9 the third of the mean field equations 
givesa ~ 91n[(l-h2,5)/(l-2^)] and n'^ ~ 49/(1-4,52)^ 
Thus, K > for sufficiently large 9 and there is a finite 
range area of phase separation. 

The thermopower S is anomalous, owing to the depen- 
dence of bandwidth on temperature through 1^(9). It 
peaks at 9c with a (5-dependent value iSniax(<5) which is 
enormous by standards of metal physics (the dimensions 
of dp/dT in Fig. ^ are restored by multiplying da/dQ 
by ks/e — 86.2 pV /K). The order of magnitude of 5, 
and the existence of a peak near Tc agree with experi- 
ments |21 . The sign of S is opposite to that predicted 
by a rigid band model (for either carrier type). At tem- 
peratures well above 9c, the temperature dependence of 
the bandwidth vanishes, and the usual sign is restored. 
Finally, at very low temperatures the thermopower also 
changes sign - the T^/^ dependence of 1 — is even- 
tually overwhelmed by the T^ dependence of p coming 
from the Sommerfeld expansion. Within our model this 
happens at a very low value of T {e.g. T/ B = 1.5 x 10^^ 
for S = 0.2). 

We have ignored so far all interactions except the dou- 
ble exchange mechanism. Electrostatic effects will in- 
hibit the full phase separation described here. The stan- 
dard RPA result X~'^{q) = X^^{q) - V{q), where Xo{q) 
is the bare charge susceptibility and V(q) = Aire'^/q'^, 
predicts an instability at finite wavevector whenever 
limq_tO Xo{q) = —V^^dx/dp is positive (V is the unit cell 
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FIG. 2. Solution to the model (|8|j9|) using an elliptical den- 
sity of states and S = | at 5 = 0.33. The bottom panel shows 
thermopower (dashed; units of fcs) and inverse compressibility 
(solid), which is negative over a band of temperatures. 



volume), with 



V 



(10) 



At low dopings, ~ a;, while at higher fillings, is pro- 
portional to the Thomas- Fermi wavevector g^F- Thus, we 
expect that electrostatic effects will give rise to a domain 
structure, at a length scale of order q~^. 

In the CMR manganites, the eg orbitals are doubly 
degenerate in the absence of JT distortions. While we 
expect intraatomic Coulomb effects to preclude multiple 
occupancy of the Cg levels, orbital ordering could signifi- 
cantly alter the picture discussed here (for example, the 
sign of the thermopower could reverse). We have also 
ignored interactions with the lattice. Our results suggest 
that the tendency towards phase separation gives rise 
to domains with a very low concentration of holes. In 
these domains, JT deformations will be favored ||2^. In 
general, the enhancement of charge fluctuations near Tc 
should induce significant lattice deformations, because 
of the different sizes of the Mn^+ and Mn^+ ions. A 
structure with domains of different electronic density and 
bandwidth enhances the tendency towards localization. 
Standard models for the infiuence of critical spin fluctu- 
ations are not applicable if, in addition, there are strong 
charge fluctuations. A domain with extra charge and 
stronger ferromagnetic order is akin to the spin polarons 
proposed in the literature. In our case, however, such a 
structure arises from magnetic effects alone. 
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